UNCLASSIFIED 


_ AD  NUMBER _ 

AD882199 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
only;  Test  and  Evaluation;  APR  1971.  Other 
requests  shall  be  referred  to  Arnold 
Engineering  Development  Center,  Arnold  ADB,  TN. 


_ AUTHORITY 

AEDC  Itr  30  Aug  1984 


THIS  PAGE  IS  UNCLASSIFIED 


AEDC-TR-71-68 


2S  197) 


RECOVERY  CHARACTERISTICS  OF  A  SINGLE-SHIELOED 
SELF-ASPIRATING  THERMOCOUPLE  PROBE  AT  LOW 
PRESSURE  LEVELS  AND  SUBSONIC  SPEEDS 


for  poBHc  release;  dfstiibntlon 

-K'  :^ccrc<L 

'O  /^uo  /9  8^ 

{Vo  mn.  ^T/'« 


C.  E.  Willbanks 
ARO,  Inc. 


April  1971 


K8ll^Bt4oaJiniited,J<J^.  S.  GovernmepoaenclgSonlyu 
covcrr8i&1e^^;^yaluationot.wltfSJ>^ardware^,/(l5m 
1971;  othpr^questr^>4ll^MCument  mus^e^^erred 
to  ApAfC  (tfJ^f^ri^‘KtterSlm-Aii^r6e  Base,  Ohio 

Arnold  ,^rt)ree.„StatiM, 


ENG/Nf  TE$t  FACILITY 

ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE 


r 


mim 


When  U.  5.  Government  drawings  specifications,  or  otlier  data  are  used  (or  any  purpose  other  than  a 
definitely  related  Government  procurement  operation,  the  Government  thereby  incurs  no  responsibility 
nor  any  obligation  whatsoever,  and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication 
or  otherwise,  or  in  any  manner  licensing  the  holder  or  any  other  person  or  corporation,  or  conveying 
any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 


Qualified  users  may  obtain  copies  of  this  report  from  the  Defense  Documentation  Center. 


References  to  named  commercial  products  in  this  report  are  not  to  be  considered  in  any  sense  as  an 
endorsement  of  the  product  by  the  United  States  Air  Force  or  the  Government. 


AEDC-TR-71-68 


RECOVERY  CHARACTERISTICS  OF  A  SINGLE -SHIELDED 
SELF -ASPIRATING  THERMOCOUPLE  PROBE  AT  LOW 
PRESSURE  LEVELS  AND  SUBSONIC  SPEEDS 


/ 


C.  E.  Willbanks 
ARO,  Inc. 


Distribution  limited  to  U.  S.  Government  agencies  only; 
covers  the  test  and  evaluation  of  military  hardware;  April 
1971 ;  other  requests  for  this  document  nrust  be  referred 
to  AFAPL  (APTP),  Wri^t-Patterson  Air  Force  Base,  Ohio 
45433,  or  AEDC  (XON),  Arnold  Air  Force  Station, 
Tennessee  37389. 


AEDC-TR-71-68 


FOREWORD 


The  ^vork  reported  herein  was  sponsored  by  the  Air  Force  Aero¬ 
space  Propulsion  Laboratory,  Air  Force  Systems  Command  (AFSC), 
Wright -Patterson  Air  Force  Base,  Ohio,  under  Program  Element 
65401F. 

The  results  of  the  research  presented  were  obtained  by  ARO,  Inc- 
(a  subsidiary  of  Sverdrup  &  Parcel  and  Associates,  Inc. ),  contract 
operator  of  AEDC,  AFSC,.  under  Contract  F40600-71-C-0002.  The 
research  was  conducted  from  September  1968  to  June  1970  under  ARO 
Project  Numbers  RD0953  and  RD0850,  and  the  manuscript  was  submitted 
for  publication  on  February  4,  1971. 

This  technical  report  has  been  reviewed  and  is  approved. 

Walter  C.  Knapp 
Lt  Colonel,  USAF 
AF  Representative,  ETF 
Directorate  of  Test 


Joseph  R.  Henry 
Colonel,  USAF 
Director  of  Test 


11 


ABSTRACT 


An  apparatus  which  produces  a  well-defined  flow  field  for  calibra¬ 
tion  of  thermocouple  probes  over  a  wide  range  of  subsonic  flow  condi¬ 
tions  at  low  pressure  levels  is  described.  Results  for  calibration  of  a 
typical  thermocouple  probe  used  in  engine  testing  are  presented.  Tem¬ 
perature  recovery  characteristics  are  presented  for  a  thermocouple 
probe  of  the  single-shielded,  self- as  pirating  type  for  total  pressures 
of  0.  75,  1.  0,  3.  0,  7.  5,  and  15  psia  with  subsonic  flow  speeds  ranging 
from’  Mach  number  0.  25  to  Mach  number  1.  0.  The  total  temperature 
of  the  air  was  approximately  670“R  for  all  the  tests.  An  analysis  of 
the  probe  recovery  characteristics  based  on  a  semi -empirical  correla¬ 
tion  was  made,  and  the  correlation  appears  to  give  satisfactory  results 
over  the  entire  Mach  number  and  total  pressure  range  considered. 
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SECTION  I 
INTRODUCTION 


The  temperature  indicated  by  a  thermocouple  probe  is  governed  by  a 
balance  between  the  net  rate  of  heat  transfer  and  rate  of  energy  storage 
at  the  thermocouple  junction.  Heat  transfer  may  occur  by  the  convection, 
radiation,  and  conduction  modes.  The  relative  importance  of  these  modes 
of  heat  transfer  depends  in  large  measure  on  probe  design.  In  the  ab¬ 
sence  of  heat  exchange  between  a  probe  and  its  environment,  a  thermo¬ 
couple  probe  usually  indicates  a  temperature  which  lies  between  the  static 
temperature  and  total  temperature  of  the  stream  in  which  it  is  immersed. 
This  indicated  temperature  is  called  the  adiabatic  probe  temperature  or 
recovery  temperature  of  the  probe.  The  difference  between  static  and 

total  temperature  is  ^  -  m2)  T.  Thus,  it  is  necessary  to  know  the 

relationship  of  the  recovery  temperature  to  static  temperature  if  the 

error  in  temperature  is  to  be  kept  below  -  -  m2)t.  Radiation  and 

conduction  heat  transfer  to  a  probe  is  not  always  negligible,  and  fre¬ 
quently,  corrections  must  be  applied  to  the  probe  indicated  temperature 
before  the  recovery  temperature,  and  finally  the  true  gas  temperature, 
can  be  obtained. 

The  relationship  of  recovery  temperature  to  total  temperature  for  a 
probe  depends  on  probe  geometry,  the  material  from  which  the  probe  is 
constructed,  probe  orientation  relative  to  the  flow  direction,  type  of  fluid, 
the  fluid  velocity,  and  the  thermodynamic  state  of  the  fluid  itself.  In 
almost  all  cases,  this  relationship  must  be  determined  experimentally. 
The  experimental  approach  is  to  immerse  a  particular  probe  in  a  well- 
defined  flow  field  and  to  compare  the  recovery  temperature  of  the  probe 
with  the  total  temperature  of  the  fluid.  By  exposing  the  probe  to  a  suf¬ 
ficiently  wide  range  of  flow  conditions,  its  recovery  characteristics  can 
be  determined  as  a  function  of  the  velocity  and  thermodynamic  state  of 
the  fluid. 

One  of  the  problems  associated  with  this  experimental  approach  is 
the  accurate  determination  of  the  total  temperature  of  the  fluid  in  .the 
calibration  apparatus.  Usually  the  total  temperature  is  measured 
directly  with  a  thermocouple  in  a  region  of  the  flow  field  where  the  fluid 
velocity  is  low.  This  technique  becomes  increasingly  less  reliable  as 
the  gas  density  is  reduced  since  stem  conduction  and  radiation  heat 
transfer  effects  become  of  greater  importance  relative  to  the  convective 
heat  transfer  in  establishing  the  thermocouple  temperature.  The 
apparatus  used  in  this  study  circumvents  the  difficult  problem  of  direct 
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measurement  of  the  fluid  temperature  in  the  apparatus.  In  the  apparatus, 
air  is  heated  at  low  velocity  in  an  indirect  heat  exchanger  having  a  large 
amount  of  surface  area  and  employing  saturated  steam  as  an  isothermal 
heat  source.  After  being  heated  to  very  near  the  steam  temperature, 
the  air  is  accelerated  in  a  nozzle  to  the  desired  Mach  number.  The 
total  temperature  of  the  fluid  stream  is,  therefore,  nearly  equal  to  the 
temperature  of  the  heat  exchanger  core  which  can  be  readily  measured. 
The  apparatus  is  well  suited  for  calibration  of  probes  at  low  gas  densi¬ 
ties  corresponding  to  high  altitude  flight  conditions. 

Figures  la  and  b  (Appendix  I)  show  a  sketch  and  a  photograph  of  a 
typical  thermocouple  probe  configuration  for  use  in  turbojet  engine  test¬ 
ing.  This  probe  is  self  aspirating  since  the  flow  around  the  probe  shield 
induces  flow  through  the  shield  and  around  the  thermocouple  junction.  In 
the  present  study,  recovery  temperature  characteristics  of  three  probes 
of  this  configuration  (Fig.  la)  were  determined  over  a  wide  range  of  sub¬ 
sonic  flow  conditions.  Saturated  steam  at  atmospheric  pressure  was 
used  as  a  heat  source;  therefore,  all  of  the  tests  were  conducted  at 
approximately  210‘’F  total  temperature. 

No  attempt  is  made  here  to  present  a  comprehensive  survey  of  the 
literature  on  thermocouple  probes.  Rather,  the  survey  is  confined 
primarily  to  those  studies  considering  probe  recovery  temperature  char¬ 
acteristics  at  low  pressures  in  subsonic  flow,  i.  e.,  characteristics  at 
low  Reynolds  numbers,  and  those  studies  on  thermocouple  probes  of  the 
self-as pirating  type  since  these  are  the  areas  of  primary  interest  in  the 
present  investigation. 

Hottel  and  Kalitinsky  (Ref.  1)  presented  the  results  of  an  investiga¬ 
tion  of  temperature  measurements  in  high  velocity  air  streams.  Three 
types  of  self-aspirating  probes  as  well  as  bare  wire  thermocouples  were 
investigated.  The  parameters  affecting  thermocouple  behavior  in  a  gas 
stream  were  noted.  Importance  of  probe  vent  area  and  thermocouple 
location  for  the  so-called  Franz  probe  was  shown,  and  temperature  re¬ 
covery  and  radiation  characteristics  were  determined  for  two  other  self¬ 
aspirating  types  of  thermocouple  probes. 

The  results  of  a  theoretical  and  experimental  study  of  thermocouple 
probe  recovery  temperature  characteristics  is  presented  in  Ref.  2  by 
Malmquist.  In  this  study,  it  was  found  that  there  exists  an  optimum  vent 
hole  area-to-internal  shield  cross-sectional  area  ratio  for  maximum 
probe  temperature  recovery,  and  the  importance  of  locating  the  thermo¬ 
couple  junction  ahead  of  the  vent  holes  was  shown.  Malmquist  also  deter¬ 
mined  that  the  material  from  which  a  probe  is  constructed  has  an  influ¬ 
ence  on  its  recovery  characteristics.  The  experiments  were  conducted 
at  a  total  temperature  near  80“ F  and  a  static  pressure  of  one  atmosphere. 
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Considerable  research  effort  was  devoted  to  the  investigation  and 
development  of  thermocouple  probes  by  NACA.  Two  reports  on  this 
research  are  particularly  pertinent  to  the  present  study.  Stickney 
(Ref.  3)  reported  time  constants  and  temperature  recovery  character¬ 
istics  for  six  thermocouple  probes.  Data  were  presented  for  a  total  tem¬ 
perature  of  approximately  80®F  for  both  subsonic  and  supersonic  Mach 
numbers  and  over  the  total  pressure  range  from  0.  2  to  2.  3  atm.  Glawe, 
Simmons,  and  Stickney  (Ref.  4)  presented  results  of  an  experimental 
study  of  thermal  radiation  and  recovery  characteristics  as  well  as  time 
constants  for  nine  probes.  Total  pressures  of  0.  2,  0.  5,  and  1.  0  atm, 
Mach  numbers  of  0.  3,  0.  6,  and  0.  9,  and  a  total  temperature  of  approxi¬ 
mately  80“F  were  the  test  conditions  in  the  determination  of  recovery 
temperature  characteristics.  One  of  the  probes  tested  is  similar  to  the 
probe  investigated  in  this  study.  In  order  that  these  two  probes  may  be 
compared  and  contrasted,  they  are  shown  in  Fig.  1. 

In  all  of  the  previously  mentioned  studies,  the  stagnation  tempera¬ 
ture  of  the  gas  stream  was  determined  through  direct  measurement  in 
each  calibration  apparatus. 


SECTION  II 
APPARATUS 


The  arrangement  of  the  experimental  apparatus  used  in  the  investi¬ 
gation  is  shown  in  Fig.  2.  It  is  connected  to  the  Engine  Test  Facility 
(ETF)  air  supply  and  exhaust  plant.  This  apparatus  is  capable  of  pro¬ 
viding  air  within  a  fraction  of  a  degree  of  the  saturation  temperature  of 
steam  at  ambient  pressure.  In  the  test  section,  velocities  up  to  sonic 
velocity  and  total  pressures  from  0.  75  to  15  psia  can  be  provided. 


2.1  TEST  CELL 

The  heart  of  the  apparatus  is  the  heat  exchanger,  a  photograph  of 
which  is  shown  in  Fig.  3.  Appendix  III  presents  an  analysis  of  the  per¬ 
formance  of  the  heat  exchanger  which  gives  the  temperature  rise  of  the 
air  through  the  exchanger  as  a  function  of  flow  rate.  After  flowing 
through  the  heat  exchanger,  the  air  enters  a  converging  nozzle  where  it 
is  accelerated  to  the  desired  Mach  number.  The  flow  from  the  nozzle 
forms  a  free  jet,  and  the  test  section  is  taken  to  be  one  inch  downstream 
of  the  nozzle  lip  which  is  well  within  the  potential  core  of  the  jet.  The 
Mach  number  in  the  test  section  is  calculated  from  the  isentropic  flow 
relations  and  the  pressure  ratio  across  the  nozzle.  The  pressure  ratio 
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across  the  nozzle  was  determined  from  the  plenum  pressure  and  the  static 
pressure  at  the  w^all  of  the  exhaust  duct  at  a  location  coincident  with  the 
nozzle  exit  plane.  This  implies  the  assumption  of  uniform  static  pres¬ 
sure  across  the  exhaust  duct  in  the  neighborhood  of  the  nozzle  exit  plane. 

As  indicated  in  Fig.  2,  the  saturated  steam  is  generated  by  boiling 
water  in  the  bottom  of  the  steam  jacket.  Heat  to  boil  the  water  is  ob¬ 
tained  by  passing  steam  from  the  steam  plant  through  copper  coils 
immersed  in  the  water.  The  saturated  vapor  is  forced  to  flow  up 
vertically  through  the  tubes  in  the  heat  exchanger  and  along  the  steam 
jacket  where  it  surrounds  the  exhaust  duct  and  nozzle  flange.  It  exhausts 
into  the  laboratory  at  ambient  pressure.  The  passages  for  flow  of  steam 
are  large,  and  flow  velocities  are  low;  therefore,  the  pressure  inside 
the  steam  jacket  of  the  apparatus  is  virtually  uniform  at  ambient  pres¬ 
sure.  Since  the  steam  is  saturated,,  the  pressure  determines  its  tem¬ 
perature;  hence,  the  steam  temperature  is  uniform  inside  the  apparatus. 

The  nominal  exit  diameter  of  the  nozzle  is  four  inches.  The  plenum- 
to-nozzle  area  ratio  is  11.  4,  which  gives  a  maximum  Mach  number  of  0.  05 
in  the  plenum.  Thus,  the  plenum  can  be  considered  to  be  a  stagnation 
chamber,  and  the  total  temperature  can  be  taken  to  be  equal  to  the  tem¬ 
perature  of  the  air  at  the  exit  of  the  heat  exchanger. 

Pressures  in  the  apparatus  were  measured  with  a  combination  of 
gages  and  differential  manometers  apropos  the  test  conditions. 

The  temperatures  were  measured  with  copper-constantan  thermo¬ 
couples.  The  electrical  circuit  for' the  thermocouple  is  shown  schemat¬ 
ically  in  Fig.  4.  The  electromotive  force  (emf)  output  for  the  thermo¬ 
couples  was  measured  with  a  potentiometer  in  a  bridge  circuit  utilizing 
an  external  electronic  galvanometer  to  indicate  a  null  condition  for  the 
circuit.  From  Fig.  4,  it  can  be  noted  that  the  thermocouples  in  the 
probes  are  referenced  to  a  thermocouple  junction  mounted  on  the  heat 
exchanger. 


2.2  TEST  ARTICLE 

A  sketch  of  the  thermocouple  probe  which  was  tested  is  shown  in 
Fig.  la.  A  photograph  of  the  probe  is  presented  in  Fig.  lb.  Three 
probes  were  tested  simultaneously  by  installing  the  probes  in  a  rake. 
The  rake  was  mounted  so  that  the  center  probe  was  on  the  centerline  of 
the  nozzle.  The  other  two  probes  were  positioned  0.  75  in.  above  and 
below  the  centerline.  All  of  the  probes  were  located  approximately  one 
inch  downstream  of  the  nozzle  lip.  The  probes  were  fabricated  from 
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stainless  steel  tubing  (Types  303  and  304),  and  the  thermocouples  in  the 
probes  were  copper -constantan.  Equality  of  aspiration  rates  for  the  three 
probes,  at  least  for  high  Reynolds  number,  was  verified  by  comparing  the 
rates  of  flow  through  the  shields  in  a  direct-connect  test.  One  end  of  a 
piece  of  plastic  tubing,  whose  other  end  was  connected  to  a  nitrogen  gas 
bottle,  was  fitted  to  the  front  end  of  a  probe  shield  with  the  vent  holes  in 
the  probe  being  left  open  to  the  atmosphere.  Thus,  nitrogen  could  flow 
from  the  bottle  through  the  tubing  and  through  the  shield  and  finally  be 
exhausted  to  the  atmosphere  through  the  vent  holes.  The  bottle  was  pres¬ 
surized  with  nitrogen,  and  the  time  required  for  the  pressure  in  the  bottle 
to  drop  a  given  amount  was  noted.  This  was  repeated  for  each  of  the 
three  probes.  The  times  for  the  three  probes  were  found  to  agree  satis¬ 
factorily,  thus  indicating  equality  of  aspiration  rates  at  high  Reynolds 
number. 

It  is  worthy  to  note  that  considerable  difficulty  was  experienced  with 
foreign  particles  becoming  lodged  in  the  vent  holes  of  the  probes  during 
testing.  Usually  there  was  a  discernible  decrease  in  the  recovery  tern- 
perature  for  the  probe  when  this  happened.  It  was  necessary  to  remove 
the  rake  from  the  cell  and  clean  the  probes  when  such  a  change  was  noted. 

Early  in  the  study,  difficulties  were  also  experienced  with  stem 
leakage  because  of  the  large  pressure  difference  across  the  stem.  This 
problem  was  solved  by  redesigning  the  probe  rake. 


SECTION  III 

RESULTS  AND  DISCUSSION 


The  recovery  temperature  characteristics  of  three  probes  were 
determined  at  a  total  temperature  of  670“R  and  total  pressures  of 
0.  75,  1.0,  3.  0,  7.  5,  and  15.  0  psia  over  a  Mach  number  range  from  0.  25 
to  1.  0.  Table  I  (Appendix  II)  presents  the  results  of  the  calibration  in  the 
form  of  the  correction  factor  A  which  is  defined  as 


A  = 


The  average  recovery  correction  factor  for  the  three  probes  tested 
is  shown  graphically  as  a  function  of  Mach  number  and  pressure  level  in 
Fig.  5.  Also  shown  in  this  figure  are  recovery  correction  factor  curves 
as  predicted  by  the  semi-empirical  correlation  developed  in  Appendix  IV. 
Figure  6  shows  smooth  curves  which  fit  the  average  recovery  correction 
factor  data  given  in  Fig.  5.  The  cross-hatched  regions  represent  the 
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spread  in  data  for  all  three  probes  tested,  after  smoothing.  Most  all  of 
the  individual  data  points  for  the  three  probes  lie  within  the  cross- 
hatched  regions. 


From  Fig,  5,  it  is  evident  that  the  recovery  correction  factor  in¬ 
creases  with  Mach  number  and  decreasing  total  pressure  level.  It  was 
found  that  generally  the  recovery  correction  factor  for  Probes  2  and  3 
showed  close  agreement,  whereas  Probe  1  gave  a  slightly  lower  re¬ 
covery  correction  factor  than  did  Probes  2  and  3,  Variation  among  the 
three  probes  is  believed  to  be  due  to  differences  in  probe  construction. 

It  should  be  noted  that  the  magnitude  of  the  variation  among  probes  is 
about  the  same  as  the  variation  observed  by  Stickney  (Kef,  3)  and 
Glawe,  Simmons,  and  Stickney  (Ref,  4), 

Also  shown  in  Figs.  6c,  d,  and  e  are  the  data  obtained  by  Glawe, 
Simmons,  and  Stickney  with  the  probe  shown  in  Fig,  Ic,  Clearly  there 
is  considerable  difference  in  the  NACA  data  and  the  data  obtained  in  the 
present  study.  This  is  not  too  surprising  in  view  of  the  differences  in 
probe  design  and  construction  materials.  Apparently,  heat  conduction 
through  the  thermocouple  leads  to  the  thermocouple  junction  plays  an 
important  role  in  recovery  characteristics  of  the  thermocouple  probe, 
and  the  effective  thermal  conductivity  of  c opper -  const  antan  thermocouple 
leads  is  much  greater  than  the  effective  thermal  conductivity  of  Chromel®- 
Alumel®  thermocouple  leads.  Moreover,  the  NACA  probe  is  about  twice 
as  large  as  the  probe  of  the  present  study  and  has  a  much  different  vent 
hole  to  internal  shield  area  ratio. 


In  Appendix  IV,  a  simplified  analysis  of  the  self- aspirating  type  of 
probe  is  presented.  The  analysis  shows  that  the  temperature  recovery 
correction  can  be  correlated  in  terms  of  free-stream  Mach  number  and 
a  Nusselt  number  index  for  the  thermocouple  inside  the  shield.  The 
results  of  the  correlation  are  shown  in  Figs.  7,  8,  and  9.  To  use  the 
correlation,  the  Nusselt  number  index 


N 


“I 


is  computed  using  Fig.  7  to  estimate  the  Mach  number  of  the  flow  past 
the  thermocouple  junction.  For  air  in  the  temperature  range  from  400 

to  700“R,  Nuj  =  V  Pq^s  where  Pq  is  in  psia.  From  the  Nusselt  number 
index,  the  corresponding  value  of  is  obtained  from  Fig,  8.  The 

value  of  Aj,gf  is  determined  from  Fig.  9  at  the  free-stream  Mach  number. 
Finally,  the  recovery  correction  A  is  determined  from  the  product 
(Aref)(^/^ref)-  dashed  lines  in  Fig.  5  were  determined  from  this 
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correlation  method.  The  correlation  is  seen  to  be  reasonably  good  over 
the  entire  Mach  number  and  pressure  range  considered. 

Although  the  correlation  cannot  be  used  to  predict  the  effect  of 
changing  probe  geometry  or  size  on  recovery  temperature  correction,  it 
can  be  used  to  predict  the  effect  of  changing  the  gas  temperature  and 
composition.  This  is  discussed  in  Appendix  IV. 


SECTION  IV 

EXPERIMENTAL  ERROR  AND  ACCURACY  OF  THE  RESULTS 


In  this  section,  the  various  sources  of  experimental  errors  and  their 
effects  on  the  results  are  discussed. 


The  recovery  temperature  correction  factor  is  defined  as 


The  temperatures  (Tq  and  Tp^^)  are  not  measured  directly  but  rather  are 
inferred  from  the  measured  quantities,  Tjj  and  (Th  -  Tpa>-  Since  the 
heat  exchanger  is  always  heating  the  air,  it  follows  that  Tq  is  less'than  Tfj 
and  Tjj  is  less  than  Tg,  although  hopefully  all  three  temperatures  are  very 
close  together.  In  Appendix  III,  the  maximum  difference  between  Tq 
and  Tg  is  determined  as  a  function  of  total  pressure  and  Mach  numbers  in 
the  test  section.  From  the  analysis,  it  follows  that 

Th  -  To  <  (Ts  -  Tbj)  exp  (-28800/G) 

where  G  is  the  rate  of  mass  flow  through  the  heat  exchanger  in  pounds 
per  hour.  It  can  be  shown  that  for  isentropic  flow  through  the  nozzle 


G  = 


M 


r-  1 
2 


cy+  i)/2Cy-  1) 


For  Pq  <  3.  0  psia,  the  difference  between  the  steam  and  entering  air  tem¬ 
perature  (Tg  -  Tb^)  was  less  than  140'’F.  ^  By  preheating  the  air, 

Tg  -  Tbj  was  kept  below  40®F  for  the  data  obtained  at  7.  5  and  15  psia. 


^With  this  temperature  difference,  the  air  leaves  within  0.  05°F  of 
the  steam  temperature  when  Pq  <  3.  0  psia. 
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This  difference  is  a  fixed  or  systematic  error  which  arises  because  of 
uncertainty  in  performance  or  efficiency  of  the  heat  exchanger.  That  is, 
the  temperature  of  the  air  leaving  the  heat  exchanger  is  not  known  exactly. 
It  is  only  known  that  it  must  be  less  than  the  steam  temperature  and  must 
be  greater  than  the  temperature  predicted  by  the  somewhat  conservative 
analysis  of  Appendix  III.  In  addition  to  this  uncertainty,  there  are  errors 
associated  with  the  measurement  of  Tjj  and  Tjj  -  Tp^^.  There  are  three 
primary  sources  of -error  in  these  quantities. 

The  temperatures  were  determined  from  thermocouples  using  the 
appropriate  standard  thermocouple  tables  without  correction.  A  single 
point  calibration  of  all  the  thermocouples  between  the  ice  point  and  boil¬ 
ing  point  for  one  atmosphere  showed  that  the  temperature  difference 
given  by  the  thermocouple  was  within  1.  25°F  out  of  178®F  in  T^,  and  if 
the  error  is  assumed  to  vary  linearly  with  temperature  difference,  it 
results  in  an  error  of  0.  007  (Th  -  Tp^^)  in  Tjj  -  Tp^^.  This  error  is  con¬ 
sistent  with  the  limits  of -errbr.-given^  by  manufacturers  for  thermocouples 
fabricated  from  "off  the  shelf"  thermocouple  wire. 

The  potentiometer  used  to  measure  the  output  of  the  thermocouple 
was  calibrated  against  a  precision  voltage  source.  The  indicated  emf  of 
the  potentiometer  was  found  to  agree  with  the  actual  emf  within  2.  5  /uv 
(0.  1"F  for  the  copper  -  cons  tantan  thermocouple  system)  over  the  range 
needed  to  measure  Th  -  Tp^^.  Agreement  was  within  5.  Omv  (0.  2°F)  over 

the  range  needed  to  measure  Th  with  respect  to  an  ice  bath  at  32®F. 

In  an  attempt  to  minimize  the  error  of  observation,  three  readings 
were  taken  at  each  test  condition.  Deviation  of  individual  readings  from 
the  average  was  less  than  ±3  nw  (0. 12“F)  for  almost  all  data  points.  By 
using  the  method  of  least  squares,  the  probable  error  of  observation, 
based  on  a  probability  of  1/2  that  an  error  taken  at  random  will  be  greater 
or  less  than  the  probable  error,  was  found  to  be  less  than  1.  5  uv  or  0.06“F 
for  virtually  all  data  points. 

The  Mach  number  of  the  flow  in  the  test  section  was  determined  from 
the  isentropic  flow  relationship 

'where  the  total  pressure  (p^)  and  the  pressure  difference  (Ap)  across  the 
nozzle  were  measured.  Since  there  are  errors  in  the  measurement  of  Pq 
and  Ap,  there  is  a  corresponding  error  in  Mach  number.  The  gages  used 
in  measuring  Pq  and  Ap  were  calibrated,  and  the  maximum  error  in 
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calculating  the  Mach  number  from  the  above  equation  including  the  ob¬ 
servational  error  is  estimated  to  be  less  than  ±0.  005  Mach  number. 

This  error  in  Mach  number  can  be  converted  to  an  error  in  A  by  taking 
the  product  of  the  slope  dA/dM  with  the  error  in  M,  namely  ±0.  005.  For 
the  present  purpose,  it  is  sufficiently  accurate  to  take  A  oc  M^,  which 
leads  to  an  additional  error  in  A  of  approximately  ±0.  001  A/M. 

If  it  is  assumed  that  the  errors  are  additive,  it  follows  that 

A  =  - 

<Th), 


with  a  maximum  positive  error  of 

(9.3A  +  A/M  +  0,236)/1000 
and  a  maximum  negative  error  of 

(9.3A  +  A/M  +  0.236)/1000  +  exp 

where  the  subscript  I  denotes  the  average  indicated  values  of  Tpj  and 
Th  “  Tpg^  from  three  readings.  This  brackets  the  estimated  maximum 
error  in  A.  This  determination  of  the  maximum  negative  error  allows 
the  possibility  of  negative  values  of  A  at  certain  flow  conditions  which  is 
an  unlikely  if  not  impossible  situation  in  reality.  The  implication  of  a 
negative  value  of  A  is  that  the  recovery  temperature  is  greater  than  the 
total  temperature  of  the  gas  stream.  For  this  reason,  the  minimum  value 
of  A  is  taken  to  be  zero. 

The  error  limits  as  given  by  the  above  expression  are  shown  on  the 
respective  data  points  in  Fig.  5.  The  random  errors  which  have  been 
accounted  for  in  the  error  analysis  contribute  only  a  small  portion  to  the 
combined  error.  Thus,  the  scatter  in  the  data  presented  in  Fig.  5  sug¬ 
gests  that  there  are  other  random  errors  in  the  data  which  have  not  been 
considered. 


^  28800^ 


SECTION  V 
CONCLUSIONS 


The  recovery  temperature  characteristics  of  a  thermocouple  probe  of 
the  self-aspirating  type  frequently  used  in  turbine  engine  testing  have 
been  determined  experimentally  over  a  wide  range  of  subsonic  flow  con¬ 
ditions.  The  apparatus  used  in  this  investigation  circumvents  the  difficult 
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problem  of  direct  measurement  of  the  fluid  temperature  in  the  apparatus. 
The  fluid  temperature  must  be  known  in  order  to  specify  the  thermo¬ 
dynamic  state  of  the  fluid  and  define  a  standard  temperature  with  which 
the  indicated  temperature  of  the  probe  can  be  compared.  In  the  apparatus, 
air  at  low  velocity  is  heated  in  an  indirect  heat  exchanger  having  a  large 
amount  of  surface  area  and  employing  saturated  steam  at  atmospheric 
pressure  as  an  isothermal  heat  source.  The  total  temperature  of  the 
fluid  stream  is  very  nearly  equal  to  the  temperature  of  the  heat  exchanger 
core  which  can  be  readily  measured.  The  apparatus  is  well  suited  for 
calibration  of  probes  at  low  gas  densities  corresponding,  for  example,  to 
high  altitude  flight  conditions. 

In  addition  to  providing  recovery  temperature  correction  factor  data 
for  a  particular  probe  configuration,  the  results  of  the  investigation  indi¬ 
cate  that  the  experimental  technique  is  sound  and  can  be  used  with  confi¬ 
dence  in  calibrating  other  probe  configurations. 

A  semi-empirical* correlation  was  developed  and  appears  to  give 
reasonably  good  results  over  the  entire  range  of  Mach  number  and  total 
pressure  considered  in  the  experimental  study.  The  correlation  permits 
the  prediction  of  the  effects  of  temperature,  pressure,  Mach  number,  and 
to  some  extent,  the  gas  composition  on  temperature  recovery  correction. 
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Shield 


Notes:  Dimensions,  wire  size,  and  wire 

type  are  the  same  as  in  General 
Electric  Drawing  No.  4012303-610. 

All  dimensions  are  In  Inches. 


a.  Schematic  of  Self-Aspirating  Thermocouple  Probe 
Tested  in  the  Present  Study 

Fig.  1  Thermocouple  Probes 
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Photograph  of  Self-Aspirating  ‘ 
Tested  in  the  Present  Study 

Fig.  1  Continue 
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Section  BE  (Not  to  Scale) 


c.  Schematic  of  Self-Aspirating  Thermocouple  Probe 
Tested  by  Glawe,  Simmons,  and  Stickney  (Ref.  4) 
Fig.  1  Concluded 
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a.  Heat  Exchanger  and  Steam  Generator 
Fig.  2  Schmatics  of  Experimental  Apparatus 
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Thermocouple  Mechanically 
to  Core  of  Heat  Exchanger 


Affixed 


Ice  Bath 


Fig.  4  Schematic  Diagram  of 


Galvanometer 


Potentiometer 


Dscd,  Potentiometer  Measures 


osed.  Potentiometer  Measures  (Tjj-Tp) 
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Recovery  Correction  Factor, 


b.  Total  Pressure  1.0  psia 
Fig.  5  Continued 


Recovery  Correction  Factor, 


c.  Total  Pressure  3.0  psia 
Fig.  5  Continued 
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Recovery  Correction  Factor, 


Free-Stream  Mach  Number,  M 

d.  Total  Pressure  7.5  psia 
Fig.  5  Continued 
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Recovery  Correction  Factor, 


0  0.2  0.4  0.6  0.8  1.0 


Free-Streara  Mach  Number,  M 

e.  Total  Pressure  15.0  psia 
Fig.  5  Concluded 
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Recovery  Correction  Factor, 


0.03 


< 


Free-Stream  Mach  Number,  M 


a.  Total  Pressure  0.75  psia 

Fig.  6  Smoothed  Recovery  Factor  Data  versus  Mach  Number  and  Total  Pressure 
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Recovery  Correction  Factor, 


Free-Stream  Mach  Number,  M 

b.  Total  Pressure  1.0  psia 
Fig.  6  Continued 
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c.  Total  Pressure  3.0  psia 
Fig.  6  Continued 
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e.  Total  Pressure  15.0  psia 
Fig.  6  Concluded 


Mach  Number  in  Shield 


Free-Stream  Mach  Number,  M 


Fig.  7  Mach  Number  of  Flow  inside  Shield  of  Probe  Shown  in  Fig.  la 
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Fig.  8  Recovery  Correction  Factor  Ratio  versus  Nusselt  Number  Index 


Reference  Recovery  Correction  Factor, 


Fig.  9  Reference  Recovery  Correction  versus  Mach  Number 
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TABLE  I 

THERMOCOUPLE  PROBE  RECOVERY  CORRECTION  FACTOR  DATA 
a.  Total  Praisura  0.75  psia 


Mach 

Number, 

M 

Total 

Pressure,  p^ 
psia 

Recovery  Correction  Factor 

Probe  1 

Probe  2 

Probe  3 

Average 

m 

0. 00049 

0. 00045 

0.  00074 

0. 000561 

■IB 

0.00157 

0.00178 

0. 00164 

0.  00166 

0.00218 

0. 00269 

0.  00290 

0.  00259 

0.402 

0.  750 

0.00315 

0.  00363 

0. 00322 

0.  00333 

0.431 

0.750 

0.00340 

0.  00391 

0.  00360 

0.  00363 

0.472 

0.  738 

0. 00515 

0. 00576 

0.  00576 

0.  00556 

0.492 

0.754 

0.00467 

0.  00618 

0. 00618 

0.  00568 

0.  553 

0.  751 

0.00736 

0.  00836 

0. 00836 

0.  00803 

0.561 

0.  750 

0. 00709 

0.0080 

0. 00779 

0.  00763 

0.591 

0.  735 

0.00808 

0.  00882 

0. 00876 

0.  00855 

0.  660 

0.753 

0.01027 

0.01153 

0. 01112 

0. 01097 

0.  735 

0.  753 

0.01309 

0.01476 

0. 01399 

0.  01395 

0.  736 

0.  749 

0.01347 

0.01467 

0.01351 

0.  01388 

0.800 

0.  732 

0.01524 

0. 01572 

0. 01593 

0.  0156 

0.866 

0.  754 

0.01558 

0. 01834 

0.  01709 

0.  0170 

0.900 

0.  748 

0. 01769 

0.  02070 

0.  01904 

0.  01914 

0.  959 

0.754 

0.01785 

0.02079 

0. 01864 

0.  01910 

Note:  Total  Temperature  is  670“R. 


b.  Total  Prassure  1.0  psia 


Mach 

Number, 

M 

Total 

Pressure,  Pq 
psia 

Recovery  Correction  Factor 

Probe  1 

Probe  2 

Probe  3 

Average 

0.  254 

1.  057 

0.  000626 

0. 000955 

0.  0009701 

0. 00085 

0.  277 

0.970 

0, 000955 

0. 00149 

0.  00173 

0.  001393 

0.  305 

0.977 

0.  000597 

0. 000896 

0.  000746 

0. 000746 

0.  348 

0.985 

0.001657 

0.  00191 

0.  00194 

0. 00184 

0.403 

1.  02 

0. 00228 

0. 00333 

0.  00333 

0.  00298 

0.407 

0.  986 

0. 00279 

0. 00327 

0.  00327 

0. 00312 

0.  493 

0.  995 

0. 00518 

0. 00525 

0.  00548 

0. 0053 

0.  537 

1.01 

0.00490 

0. 00603 

0.  00603 

0.  00565 

0.  591 

0.992 

0.  00755 

0. 00766 

0.  00788 

0.  0077 

0.  650 

1.  004 

0.00937 

0. 00996 

0.  00973 

0. 00969 

0.  713 

1.  00 

0. 00825 

0. 00975 

0.  00966 

0.  00922 

0.  743 

1.01 

0.01164 

0. 01222 

0.  01203 

0.01196 

0.  770 

1.01 

0.00991 

0. 01134 

0.  00111 

0.  01080 

0. 8036 

0.998 

0.0139 

0. 0149 

0.0145 

0.  0144 

0.987 

1.01 

0.01679 

0. 01896 

0.  0183 

0.  01802 

1.  01 

1.  01 

0.01597 

0.01788 

0.01772 

0.0172 

Note:  Total  Temperature  is  670®R. 
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TABLE  I  (Concluded) 
c.  Total  Pressure  3.0  psia 


Mach 

Number, 

M 

Total 

Pressure, 

psia 

Recovery  Correction  Factor 

Probe  3 

Average 

IE9 

3.02 

0.000119 

0.  000239 

0.  000269 

0.  000209 

3.  00 

0.00228 

0.  00260 

0.  00260 

0. 00249 

■StI 

2.  997 

0.00548 

0.  00625 

0.  00603 

0. 00592 

■19 

2.  980 

0.00970 

0.  OHIO 

0.  01030 

0. 0104 

3.  00 

0.01306 

0. 01572 

0.01443 

0. 0144 

0.  904 

3.00 

_ 

0. 01439 

0.01646 

0. 01482 

0.  0152 

Note:  Total  Temperature  is  STO^R. 


d.  Total  Pressure  7.5  psia 


Mach 

Number, 

M 

Total 

Pressure,  p^ 
psia 

Recovery  Correction  Factor 

Probe  1 

Probe  2 

Probe  3 

Average 

0.  248 

7.  48 

.... 

. 

0.  292 

7.54 

0.  00025 

0.  00037 

0. 00045 

0.  000358 

0.  320 

7.  50 

0. 00051 

0.  00045 

0.  000037 

0. 000185 

0.  390 

7.497 

0.  000134 

0.  00149 

0. 00085 

0.  000826 

0.  440 

7.  507 

0.  002104 

0,  002940 

0. 002746 

0.  00260 

0.  621 

7.49 

0. 00452 

0.  00581 

0. 00588 

0.  00540 

0.  649- 

7.  49 

0.  00445 

0.  00633 

0.  00576 

0.  00551 

0.  796 

7.  50 

0. 00734 

0.  00975 

0. 00891 

0. 00867 

0.  893 

7.  50 

0.  00885 

0.  01264 

0. 01148 

_ 

0.  01099 

Note:  Total  Temperature  is  eTO^R. 


e.  Total  Pressure  15.0  psia 


Mach 

Number, 

M 

Total 

Pressure,  p^ 
psia 

Recovery  Correction  Factor 

Probe  1 

Probe  2 

Probe  3 

Average 

0.  403 

15.  0 

0.  000858 

0.  000896 

0. 00045 

0. 000715 

0.  503 

15.  0 

0. 00119 

0. 00196 

0. 00142 

0. 00152 

0.  597 

15.0 

0.00218 

0.  00309 

0. 00248 

0. 00258 

0.770 

15.0 

0. 00399 

0. 00491 

0. 00476 

0. 00455 

Note:  Total  Temperature  is  670®R. 
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The  purpose  of  this  appendix  is  to  obtain  an  estimate  of  the  differ¬ 
ence  between  the  air  leaving  the  heat  exchanger  and  the  steam  tempera¬ 
ture. 


Figure  3  shows  a  photograph  of  the  fin-tube  heat  exchanger  used  in 
the  apparatus  to  heat  the  air  to  a  known  temperature.  The  heat  exchanger 
has  an  effective  length  of  approximately  11  ft  and  a  1-ft -square  cross 
section.  The  aluminum  fins  are  0.0125  in.  thick,  and  there  are  approxi¬ 
mately  110  fins  per  foot.  The  brass  tubes  have  an  outside  diameter  of 
5/8  in.  with  a  1/32-in.  wall  thickness  and  are  arranged  in  equilateral, 
staggered  triangles  on  1-1 /2-in.  centers. 

The  rate  of  heat  transfer  to  the  air  from  the  steam  for  a  differential . 
length  of  heat  exchanger  (dx)  is 

dO  =  (Uf  Wp  +  Ut  Wt)  (Ts  -  Tb)  dx  (III-l) 

where  Up  and  Ut  are  the  overall  heat  transfer  coefficients  for  the  fins 
and  tubes,  respectively,  and  Wp  and  Wt  are  surface  areas  per  unit 
length  consistent  with  the  definitions  of  Up  and  Up.  Conservation  of 
energy  requires 

dQ  =  GCpdTb  (III-2) 


Since  saturated  steam  at  constant  pressure  is  used  as  the  heat 
source,  the  steam  temperature  T5  is  constant.  If  it  is  assumed  that 
Up,  Ut,  and  Cp  are  constant,  Eq,  (III-2)  can  be  integrated  to  give 


^  Tbjj-T,  (UfWk  +  UtWt)L 
^  Tbi-T.  '  Gc;; 


(III -3) 


where  Tb2  and  Tbj  are  inlet  and  exit  air  temperatures,  respectively, 
and  L  is  the  total  heat  exchanger  length. 


The  heat  transfer  coefficients  Up  and  Ut  must  be  evaluated  through 
consideration  of  airflow  in  the  heat  exchanger,  condensation  of  steam  in¬ 
side  the  tubes,  and  heat  conduction  through  the  fins  and  tubes.  It  can  be 
shown  that 

"  ''f  -  1  ^  D.  ^  D./D.  ^  D.  <ni-4) 

hf  .  2  k']'  hg  Dj 
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where  h-p  is  a  suitable  heat  transfer  coefficient  for  flow  over  a  finned 
tube  and  hg  is  the  heat  transfer  on  the  steam  side  of  the  tubes.  The  sur¬ 
face  area  Wx  is  the  total  surface  area  of  the  tubes  alone  exposed  to  the 
airflow  per  unit  length  of  heat  exchanger.  The  heat  transfer  coefficient 
for  the  fins  can  be  shown  to  be 


1  Ap  ^  Pp/Dj  Af  (III-5) 

hp  2JrkFt  pDit  hs 

where  Ap  is  the  surface  area  of  fin  exposed  to  the  air  per  tube  and  per 
fin.  For  this  case,  the  surface  area  Wp  is  the  total  fin  area  per  unit 
length  of  heat  exchanger  exposed  to  the  airflow. 


In  order  to  proceed  further  with  the  analysis,  it  is-necessary  to 
have  ^stimates  of  the  heat  transfer  coefficients  hg,  hp,  and  h'p  as  well 
as  the  fin  efficiency  rjp. 


From  data  on  condensation  of  saturated  steam  in  a  vertical  tube  pre¬ 
sented  in  Ref.  5,  it  can  be  inferred  that 


hs 


Btu 


(III-6) 


when  the  vapor  velocity  in  the  tube  is  low.  In  this  formula,  H  is  the 
cooled  height  of  the  tube  which  is  one  foot  for  the  heat  exchanger  under 
consideration.  Thus, 

hs  =  1065  Btu/ft2-hr-°R 


In  cross  section  along  the  direction  of  airflow,  the  fins  in  the  heat' 
exchanger  have  a  wavy  surface  to  promote  turbulence  in  the  flow  for 
effective  heat  exchange.  However,  for  purposes  of  analysis,  the  fins 
will  be  taken  to  be  flat  and  flow  between  any  two  adjacent  fins  will  be 
treated  as  laminar  and  fully  developed  channel  flow.  Clearly,  this  will 
yield  a  conservative  estimate  of  hp.  For  fully  developed  laminar  flow 
in  a  channel.  Ref.  6  gives 

2hF(S-c) 

- k -  "" 

*“a 

Since  the  thermal  conductivity  of  air  kg^  does  not  change  markedly  over 
the  temperature  range  of  interest 

hF  =  7.9  Btu/ft^-hr-Tl 
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In  view  of  the  uncertainty  in  calculating  the  heat  transfer  coefficient 
h-p  for  a  fin  tube  bank,  hp  will  be  taken  to  be  equal  to  hp.  This  is 
partially  justified  by  the  fact  that  the  total  fin  area  is  an  order  of  mag¬ 
nitude  greater  than  the  total  tube  area  exposed  to  the  air.  Thus, 

Ht  =  7.9  Btu/hr-ft2-°R 

An"  estimate  of  the  fin  efficiency  Hp  can  be  inferred  from  the  solu¬ 
tion  for  heat  conduction  in  a  circular  fin  of  rectangular  profile  attached 
to  a  circular  tube  if  the  equivalent  diameter  of  the  circular  fin  is  taken 
to  be  the  distance  between  centers  of  the  heat  exchanger  tubes.  This  is 
a  solution  given  in  many  elementary  texts  on  heat  transfer.  From  the 
numerical  results  given  in  Ref.  7,  it  can  be  estimated  that 

TfF  =  0.9 


Taking  the  conductivity  of  the  brass  tubes  to  be  35  Btu/hr-ft-‘’F  and 
substituting  the  appropriate  dimensions  and  heat  transfer  coefficients 
into  Eqs.  (III-4)  and  (III-5)  give 

Ut  =  7.84  Btu/hr-ft®-°F 


and 


Up  =  3.3  Btu/hr-ft^'-'^R 


There  are  684  tubes  in  the  heat  exchanger  and  a  total  1855  ft^  of  fin  sur 
face  area;  thus 


and 


Wt  =  9.lftVft 
Wp  =  169ftVft 


It  follows  that 


(Tbj  -  Ts)  =  (Tbi  -  Ts)  expjl 


6904  Btu/hr-°F 
GC„ 


The  left-hand  side  of  this  expression  gives  the  difference  between  the 
temperature  of  the  air  leaving  the  heat  exchanger  and  the  steam  tem¬ 
perature.  The  rate  of  mass  flow  through  the  heat  exchanger  is  governed 
by  the  total  pressure,  total  temperature,  nozzle  exit  area,  and  Mach 
number  in  the  test  section. 


I 
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APPENDIX  IV 

AN  ANALYSIS  OF  A  SELF-ASPIRATING 
THERMOCOUPLE  PROBE 


The  temperature  indicated  by  a  thermocouple  probe  at  steady  state 
is  governed  by  a  balance  between  the  conduction,  convection,  and  radia¬ 
tion  modes  of  heat  transfer  at  the  thermocouple  junction.  The  complexity 
of  the  flow  around  the  thermocouples  in  the  probes  (Fig.  1)  makes  it  im¬ 
possible  to  analyze  them  in  a  rigorous  manner  and  to  predict  their 
recovery  characteristics  theoretically.  However,  it  was  found  that  a 
reasonably  good  semi -empirical  correlation  of  the  recovery  temperature 
correction  can  be  made.  The  correlating  parameters  are  derived  below. 
Radiation  heat  transfer  to  the  thermocouple  junction  in  a  self -as pirating 
probe  of  the  type  shown  in  Fig.  1  may  be  neglected  in  comparison  with 
the  other  two  modes  of  heat  transfer.  If  there  were  no  heat  conduction 
down  the  thermocouple  leads  to  the  stem,  the  junction  temperature  (Tp^^) 

would  approach  the  total  temperature  of  the  gas  (Tq)  since  the  Mach  num¬ 
ber  in  the  shield  (Mg)  is  always  small.  Since  there  is  conduction  down 
the  leads,  however,  the  equilibrium  temperature  of  the  junction  will  lie 
between  the  total  temperature  of  the  gas  and  a  temperature,  T-pB  say, 
characterizing  the  temperature  in  the  region  where  the  thermocouple  ' 
leads  enter  the  stem. 


From  a  dimensional  analysis,  it  follows  that  (Tp^  -  Tq)/(Tpb  -  Tq) 
should  be  a  function  of  only  the  parameter  h^/k  where  h  is  the  average 
convective  heat  transfer  coefficient  for  the  thermocouple,  k  is  an  effec¬ 
tive  thermal  conductivity  of  the  thermocouple  leads,  and  is  a  character¬ 
istic  dimension  of  the  thermocouple  junction.  The  smallness  of  the 
thermocouple,  wire  and  velocity  inside  the  shield  yields  a  low  Reynolds 
number  and  laminar  flow  over  the  thermocouple.  It  is  well  known  that 
the  Nusselt  modulus  for  heat  transfer  in  laminar  flow  varies  as  the 
square  root  of  the  Reynolds  number  and  the  cube  root  of  the  Prandtl  num¬ 
ber  for  body  shapes  ranging  from  the  flat  plate  to  the  sphere.  It  is, 
therefore,  not  unreasonable  to  assume  that  this  dependence  holds  for  the 
Nusselt  modulus  of  the  thermocouple.  Thus, 


where  a  is  a  constant  of  proportionality  and  the  gas  properties  are  eval¬ 
uated  at  the  total  conditions  of  the  free  stream. 
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An  estimate  of  the  Mach  number  inside  the  shield  can  be  obtained  by 
assuming  that  the  flow  expands  isentropically  to  the  free -stream  static 
pressure  at  the  vent  holes.  It  follows  that 

where  Cj^  is  the  discharge  coefficient  for  the  vent  holes.  The  discharge 
coefficient  for  a  sharp-edged  orifice  is  approximately  constant  with  a 
value  of  0.  6  for  a  Reynolds  number  greater  than  15  based  on  orifice 
diameter.  Although  the  flow  situation  for  the  vent  holes  is  somewhat 
different  than  the  usual  orifice  flow  situation,  it  is  probably  sufficiently 
accurate  to  take  Cp  =  0.  6  in  Eq.  (IV-2).  The  vent  hole -to -shield  ratio 
for  the  probe  shown  in  Fig.  la  is  approximately  0.  183.  The  internal 
shield  Mach  number  (Mg)  is  plotted  versus  M  in  Fig.  7.  Although  k  and  i 
are  unknown,  they  are  fixed  for  a  given  probe.  It  is  convenient  to  com¬ 
bine  a,  k,  and  i  into  a  single  dimensional  constant  Because  the 

exact  value  of  this  constant  is  imimportant  in  the  correlation,  it  was 
chosen  to  make  hi/k  equal  to  unity  when  Pq  =  1  atm,  Tq  =  670“R,  and 
M  =  1.0.  From  Fig.  7,  it  can  be  determined  that  Mg  =  0.  0635  when 
M  =  1.  0.  Hence,  ay[T/k  =  1.  1  x  10^  ft^/^-hr-^R/Btu. 

Since  the  stem  temperature  (Tg^)  is  unknown,  (Tp^^  -  Tq)/(Tbx  -  Tq) 
cannot  be  correlated  directly  with  Hi /Ic.  Both  experiment  and  theory  indi¬ 
cate  that  the  recovery  temperature  for  a  wide  range  of  geometrical  shapes 
does  not  depend  on  Reynolds  number  except  during  the  transition  from 
laminar  to  turbulent  flow.  Therefore,  it  is  not  without  justification  to 
assume  that  Tgx  is  independent  of  Reynolds  number.  It  should  be  re¬ 
marked  that  the  stem  length  was  many  stem  diameters  long  and  that  a 
calculation  based  on  theory  of  heat  conduction  in  extended  surfaces  showed 
stem  conduction  to  be  negligible.  It  follows  that,  at  a  given  Mach  number 
and  total  pressure, 

A  ftiif/k] 

where  the  subscript  ref  refers  to  values  of  A  and  hi  Ik  at  the  given  Mach 
number  and  a  reference  total  pressure. 

Figure  8  shows  A/Aj.gf  versus  hi  Ik  for  the  data  given  in  Table.  I  for 
a  reference  total  pressure  of  1  psia.  The  values  of  Aj,gf  were  obtained 
from  the  curve  given  in  Fig.  9,  which'was  obtained  by  drawing  a  smooth 
curve  through  the  average  recovery  correction  for  the  three  probes  at 
1  psia.  The  smooth  curve  in  Fig.  8  was  taken  to  represent  the  data 
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plotted  in  Fig.  IV- 1.  From  Fig.  8,  it  can  be  noted  that  A/^ref  shows  a 
weak  dependence  on  h  i  /Ic  at  small  values  of  h  i /It.  Thus  to  a  good 
approximation  A/A^q^  can  be  correlated  in  terms  of  hi/Ic  alone;  that  is, 

the  variation  with  (hj2/k)i.ef  be  neglected. 

In  order  to  obtain  A  from  the  correlation,  ^ref  is  obtained  from 
Fig.  9  and  A/Aj-ef  is  obtained  from  the  smooth  curve  in  Fig.  8.  The 
product  of  Aj.gf  and  A/Aj-gf  is,  of  course,  A.  The  dashed  curves  shown  in 
Fig.  5  were  obtained  in  this  way.  It  can  be  seen  that  the  correlation  is 
reasonably  good  over  the  entire  range  of  total  pressure  and  Mach  num¬ 
ber  considered. 

It  is  instructive  to  consider  the  effect  of  total  temperature  on  the 
recovery  correction  factor  as  predicted  by  the  correlation.  By  assuming 
that  the  viscosity  and  the  thermal  conductivity  of  air  vary  as  the 
0.  7  power  of  the  absolute  temperature,  it  follows  that 

k 

/ 

Clearly,  changing  Tq  by  any  reasonable  amount  would  influence  A  very 
little.  It  would  be  desirable  to  conduct  experiments  over  a  range  of  tem¬ 
perature  to  confirm  this  conclusion. 

In  order  to  obtain  further  confidence  in  the  correlation  technique, 
the  data  for  the  self-aspirating  probe  {Fig.  Ic)  presented  in  Ref.  4  were 
correlated,  and  the  results  of  the  correlation  are  presented  in  Fig.  IV-2. 
The  correlation  is  good  for  the  NACA  probe.  It  should  be  remarked  that 
the  correlation  curves  given  in  Figs.  7,  8,  and  9  are  not  applicable  to 
the  NACA  probe  because  the  curves  were  obtained  specifically  from  the 
recovery  correction  data  on  the  probe  of  the  present  study.  The  cor¬ 
responding  correlation  curves  for  the  NACA  probe  were  found  to  differ 
considerably  from  the  curves  in  Figs.  7,  8,  and  9. 
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Recovery  Correction  Factor, 


Fig.  IV-2  Recovery  Correction  Factor  for  NACA  Probe  (Probe  6,  Ref.  4) 
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was  made,  and  the  correlation  appears  to  give .satisfactory  results  over 
the  entire  Mach  number  and  total  pressure  range  considered. 
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